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Abstract 

In this four-part prospectus, we first give a brief introduction to the motivation for studying 
entanglement entropy and some recent development. Then follows a summary of our recent work 
about entanglement entropy in states with traditional long-range order. After that we demonstrate 
calculation of entanglement entropy in both one-dimensional spin-less fermionic systems as well as 
bosonic systems via different approaches, and connect them using one-dimensional bosonization. 
In the last part, we briefly sketch the idea of bosonization in high-dimensions, and discuss the 
possibility and advantage of approaching the scaling behavior of entanglement entropy of fermions 
in arbitrary dimensions via bosonization. 
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1. INTRODUCTION 



Entanglement is the hallmark as well as most counter-intuitive feature of quantum me- 
chanics. Its contradiction with one of the most important concepts in classic physics - locality 
- has caused skepticism and controversy[14] ever since quantum mechanics was formulated. 
It was only after the genuine insight of Dr. John Bell [5], where based on several general 
assumptions he derived a set of inequalities for two physical observables that any local the- 
ory should obeythat the difference between truly quantum entanglement and classic or local 
correlation can be tested and distinguished experimentally. However, the experiments only 
become viable after decades[17] since these inequalities were formulated. The overwhelming 
majority of experiments done nowadays support quantum entanglement [30]. Even though 
there are still critics pointing out that these Bell test experiments are not problem-free, 
the existence of such experimental problems which are often referred to as 'loopholes' may 
affect the validity of the experimental findings, they are still demonstrating that quantum 
entanglement is physical reality. 

Entanglement has attracted more attention due to the development of quantum informa- 
tion and quantum computation science[29] where entanglement is considered major resource 
of quantum information. 

Among various ways of quantifying entanglement, bipartite block entanglement entropy 
has emerged as a concept of central importance, not only in quantum information science, 
but recently in other branches of physics as well. In condensed matter/many-body physics 
the entanglement entropy has been increasingly used as a very useful and in some cases indis- 
pensable way to characterize phases and phase transitions, especially in strongly correlated 
phases [1]. 

This bipartite block entanglement entropy is defined as the following. For any given 
density operator of a pure state p, we can divide the system into two parts, then we partially 
trace out either part from the division: 

Pi = tr 2 (p), p 2 = tri(p). (1) 

And the entanglement entropy is defined as von Neumann entropy of the reduced density 
matrix: 
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E = Ei = Ei 



tr(pi lnpi) 



tr(p 2 lnp 2 ) 



(2) 



In this context the most important result is perhaps the so-called area law [6, 35], which 
states that in the thermodynamic limit, the entropy should be proportional to the area of 
the boundary that divides the system into two blocks. There are a very few important 
examples [1, 19, 40] in which the area law is violated, most of which involves quantum 
criticality[7, 9, 15, 21, 33, 34, 36]; the specific manner with which the area law is violated 
is tied to certain universal properties of the phase or critical point. In some other cases, 
important information about the phase can be revealed by studying the leading correction to 
the area law [16, 25, 27]; for example this is the case for topologically ordered phases[25, 27]. 

It is to our particular interest to study the entanglement entropy of many-fermion systems 
which are most encountered in condensed matter physics. There have been numerous works 
on this topic, mostly numerical. There are also some analytic results in ID spin-less fermions, 
based on conformal field theory argument [9, 21] or fermionic wave functions[23, 24], giving 
entanglement entropy | logL where L is the subsystem size. For higher dimensions, Wolf [40] 
establishes a relation between structure of the Fermi sea and the scaling of the entropy for 
a finite nonzero Fermi surface: 



while Gioev and Klich[19] provide an explicit geometric formula for the entropy of free 
fermions in any dimension d as L — > oo by making a connection with Widom's conjecture [39]: 



1.1. General Form of Reduced Density Matrix 

In this section, we will consider a general formalism for obtaining the reduced density 
matrix in Fock space. 

Consider a general pure state \^q), and we want to divide it into two parts in some 
orthonormal basis. Let the basis be \B), and the basis of the two divisions are denoted as 
|e), \s) respectively. Here we use e indicating the environment which we want to trace out 
and s denoting the subsystem of interest. Then in general we can write: 



S ~ L logL. 



(3) 
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I^o) = ^c B \B) = ^c e , s |e) ® \s). (5) 

B e,s 

So the density matrix of the whole system is written as: 



Po = l^oX^ol = I Yl ® I ( <2,s 2 {^\ ® (s 2 \ 1 ■ (6) 

\ei,si / \e2,S2 / 

In order to easily trace out the environment part, we collect the states with the same 
environment parts, which yields: 



Po 



v ei si / \ e2 S2 J 



Then we move on to the calculation of the reduced density matrix by tracing out the 
environment: 



Ps = tr e (p ) = <eo|(X)(I>i,.il a i>) ® |ei»0BE c U^|) S (e 2 |)|e ) • (8) 

eo \ ei si e2 S2 / 

Since |e) is chosen to be orthonormal, so (e |ei) = S e0j6l . the above equation is reduced 



to: 



p s = tr e (p ) = ( (J] Ce o^l Sl ))E c eo^( s 2|) J 



(9) 



sx,S2 eo 



So from Eq.9 we can get that the matrix element of the reduced density matrix of sub- 
system s is: 

(Ps)( Si ,s 2 ) = (Sl|p s |s 2 ) = Y C ^i C ls 2 ( 10 ) 

e 

According to Eq.5, we can write: 

(ps)( Sl ,s2) = ^c eiS1 c* S2 = (*o|aJ 2 a«i|*o), (11) 
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where a s , which we shall call 'state annihilation operator', is defined within the Hilbert 
space of the subsystem as a Sa \s) = 6 SQ>S \0), |0) being the vacuum state of the subsystem. 
The specific form of a so need be addressed in specific systems. 

Now we can write the reduced density matrix in a compact form: 



This formalism is general and works for any complete basis of the Fock space of concern. 

As long as we are only concerning about the density matrix of the subsystem, we can 
replace our 'state annihilation operators' by any complete set of operators that span the 
complete Hilbert space of the subsystem. 

2. ENTANGLEMENT ENTROPY IN STATES WITH TRADITIONAL LONG- 
RANGE MAGNETIC ORDER 

A lot of works have been done on entanglement entropy in various aspects, however, 
there have been relatively few studies of the behavior of entanglement entropy in states 
with traditional long-range order[3, 26, 37]. This is perhaps because of the expectation 
that ordered states can be well described by mean-field theory, and in mean-field theory the 
states reduce to simple product states that have no entanglement. In particular in the limit 
of perfect long-range order the mean-field theory becomes "exact", and the entanglement 
entropy should vanish. In this section, we will summarize our recent work[13] on those states 
and show that this is not the case, and interesting entanglement exists in states with perfect 
long-range order. We will study two exactly solvable spin-1/2 models: (i) An unfrustrated 
antiferromagnet with infinite range (or constant) antiferromagnetic (AFM) interaction be- 
tween spins in opposite sub-lattices, and ferromagnetic (FM) interaction between spins in 
the same sub-lattice; (ii) An ordinary spin-1/2 ferromagnet with arbitrary FM interaction 
among the spins[32] 1 . While the ground states have perfect long-range order for both mod- 
els, we show that they both have non-zero entanglement entropy that grow logarithmically 

1 The entanglement property of this model has bee previously studied in Ref. [32]. Here we introduce 
a different definition of the entanglement entropy that properly takes into account the ground state 
degeneracy, and can be calculated much more straightforwardly using symmetry consideration. See Sec. 




(12) 



2.3. 
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FIG. 1: Two-sub-lattice model: two sub-lattices labeled A and B interpenetrating each other, 
with the size of the subsystem. 



2.1. Antiferromagnetic Spin Model and the Ground State 

We consider a lattice model composed of two sub-lattices interpenetrating each other 
as in Fig. 1, with interaction of infinite range, i.e., every spin interacts with all the other 
spins in the system, with interaction strength independent of the distance between the spins. 
Within each sub-lattice, the interaction is ferromagnetic, and between the sub-lattices the 
interaction is antiferromagnetic; as a result there is no frustration. The Hamiltonian is 
written as, 

H = —J A ^ Si ■ Sj — J B ^ Si ■ Sj + J ^ Si ■ Sj, (13) 

i,j<EA i,j&B ieA,j£B 

with Ja-iJbiJo > 0. The ground state of (13) can be solved in the following manner [42]. 
Introduction the following notations 



Sj, Ss — '^^Si, S — Sa + Sb 



(14) 



ieA ieB 
then the Hamiltonian can be written as, 

H = —JaS 2 a — JbS b + JqSa ■ Sb 

= -JaS\ - JbS 2 b + y (£ 2 - S\ - S%) 
Then it is not hard to see that the ground state should be: 

\S A S B ; Sm) 

For simplicity we only consider the simplest case with Na = Nb = N, thus the total 
system size is 2N. Then the ground state is reduced to an antiferromagnetic ground state 
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N A 


N B 


N A 


Nb 


2 ■ 


2 ' 


2 


2 



m). (15) 



which has zero total spin, 

\S A S B] Sm) = Iff; 00). (16) 
This ground state has perfect Neel order, as manifested by the spin-spin correlation function, 

(Si ■ Sj) = i,j e A or i,j e B; 

i i i < 17 > 

(Si ■ S,) = -- - ^ -> --, i£A and j € B. 
2.2. Reduced density matrix and Entanglement Entropy 

We divide the system spatially into two subsystems which are labeled 1 and 2 respectively 
and study the ground state entanglement entropy E between these two subsystems. 

Following is a brief summary about how to solve for the explicit form of the reduced 
density matrix. 

First, we further decompose the system into four parts, S Al , Sa 2 , Sbx, Sb 2 , with 



= EieAAiei &u ^a 2 — EieAAie2 Si] 

SBi — EieBAiGl ^B 2 = EjgBAie2 

Therefore, these operators satisfy the following relations, 

Sa = Sa± + Sa 2 , Sb = Sbx + Sb 2 ] 
S\ = Sa x + Sbx, S2 = Sa 2 + Sb 2 - 



(18) 



(19) 



Here we note that, as discussed in the previous section, the spin state within each sub-lattice 
is ferromagnetic. This means that not only must the total spin quantum numbers of S A and 
Sb take their maximum values, but the total spin quantum numbers of S Al , Sa 2 , S Bl and 
Sb 2 must also take their maximum values. More importantly, these values are thus fixed, 
which enables us to treat the operators S Al , Sa 2 , S Bl and Sb 2 as four single spins, and in 
what follows we shall denote these operators by their corresponding spin quantum numbers. 
The problem is then that we are given a four spin state in which the spins Sa 1 and S 1 ^ are 
combined into a state with total spin S A and the spins S Bl and S B2 are combined in a state 
with total spin Sb, and then these two states are combined into a total singlet (resulting in 
the ground state of our long-range AFM model), and we must express this state in a basis 
in which the spins S Al and Sa 2 are combined into a state with total spin Si and Sb x and 
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Sb 2 are combined into a state with total spin S2. This change of basis involves the familiar 
LS-jj coupling scheme. At this point, to obtain the reduced density matrix we only have to 
re-express the density matrix in bases of (5 , imi5 , 2m 2 ; Sm\ by means of LSjj coupling, then 
trace out either (Simi\ or (5 , 2 m 2 |. 
After some algebra we arrive at: 

Pi = tr {2 )(p) = ^ |A 5imi | 2 |5'imi)(S'imi|. (20) 

Si mi 

where A Sl , mi is given by 

(N + 1)(N - JVi)!JVi! 



A 



Si, mi 



X 



_(N — Ni — S X )\(N - iVi + S 1 + 1)! 
(iV — iVi)!iVi! 



(21) 



(7Vi-5'i)!(^i + S , i + l)!. 
The bipartite entanglement entropy between subsystems 1 and 2 is then given by, 

E = E l = - |A Simi | 2 ln(|A 5imi | 2 ). (22) 

5i,mi 

Here we note that, although Xs imi is written with an explicit mi dependence, the actual 
expression is independent of As a result, we can eliminate the summation over mi from 
(23) by multiplying by a factor of 2S\ + 1. The final expression for the entanglement entropy 
is then, 

E = E l = -J2( 2S i + l)|Ac ?1 | 2 ln(|A. ?1 | 2 ). (23) 

Si 

Next we shall present the asymptotic behavior of this bipartite entanglement entropy in 
two limiting cases: 

• N\ — N2 — N, this case gives the saturated entropy at fixed N since intuitively E 
should increase with the subsystem size. 

• 1 Ni N, in this limit we are considering system's entanglement with its (much 
larger) environment, and generically we should be able to find that the entropy should 
be independent of the total system size as N — > 00, which is indeed what we find. 

By extracting the asymptotic behavior of the eigenvalues of the reduced density matrix 
\s 1 , turning the summation into integral and forcing the normalization condition of those 
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eigenvalues, the entanglement entropy can be analytically calculated: 

E x ~- f A(2S 1 + l) e - 5 ^iv^vr + ^ 
Jo 

C2( 1 1 



x In ^Ae~° lV7 ^^'j dS! ( 2 4) 
ln[Ni _N? , lJn(N-N l)Nl 



N 2 V JV 

First let us consider the equal partition case, which presumably gives the upper limit 
of the entanglement entropy for a given total system size. For simplicity we set N to be 
even (note that the total system size is 2N), thus N 1 = ^, and 4 = |(7V + V Nir). So the 
entanglement entropy becomes, 

N 1 



E = ln h -VTriV 

4 4 / (25) 

~ In - In 4 ~ In - 1.38629. 
For the unequal partition case, 1 C iVi < iV, we can expand the entropy as follows, 



E~\n 



N 




lniV 1 + ln|l-^ + i,/^i^l (26) 
Ni I / A 7 ! x 2 



From this expression we see that when the assumed condition is satisfied the entropy indeed 
depends only on the subsystem size to leading order. 

2.3. Ferromagnetic model and its entanglement entropy 

In this section we consider a ferromagnetic (FM) spin-1/2 model on an arbitrary lattice 
with N sites, 

H=-^J ij S i -S j , (27) 

with Jij > 0. The ground state is the fully magnetized state \SM) with S = N/2 and 
M = —S, —S + !,■•■ , S, and is clearly long-range ordered: (Si ■ Sj) = 1/4. However there is 
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a crucial difference between the FM ground state and the AFM ground state studied earlier: 
the FM ground state has a finite degeneracy, and thus the system exhibits a non-zero entropy 
even at zero temperature, E = \og(2S+ 1) = log(iV+ 1), resulting from the density matrix 
of the entire system, 

s s 
P = xf-j E \SM)(SM\ = -J— \SM)(SM\ (28) 

M=-S M=-S 

In this case the entanglement entropy between two subsystems (1 and 2) is defined in 
the following manner. We first obtain reduced density matrices for subsystems 1 and 2 by 
tracing out degrees of freedom in 2 and 1 from p: 

s 



Pi = tr {2) p = ^— E tr {2) (\SM)(SM\) 

1 5 

= w^i E 

M=-S 

1 5 

P2 = tr (1)P = w - Y E tr {l) (\SM)(SM\) 

M=-S 

1 3 

-—j E p2hu 



(29) 



N +- 

and calculate from them the entropy of the subsystems, E\ and E 2 . The entanglement 
entropy is defined as[12] 2 

E = (E 1 + E 2 -E )/2. (30) 

For the present case E\ and E 2 can be easily obtained from the following observations, (i) 
Because the total spin is fully magnetized, so are those in the subsystems: S\ = N\/2 and 
S 2 = N 2 /2. Thus this is a two-spin entanglement problem, (ii) Because the total density 
matrix p is proportional to the identity matrix in the ground state subspace, it is invariant 



2 This definition is probably not unique. It is one half of the "mutual information" introduced in Rcf. 
12, 41, and reduces to Eq. (2) when p is that of a pure state. The same definition was used by Castelnovo 
and Chamon [Claudio Castelnovo, and Claudio Chamon, Phys. Rev. B 76, 174416 (2007)]. The name 
"mutual information" may be first coined by Adami and Cerf [C. Adami and N.J. Cerf, Phys. Rev. A 
56, 3470 (1997)] and Vedral, Plenio, Rippin and Knight [V. Vcdral, M.B. Plenio, M.A. Rippin, and P.L. 
Knight, Phys. Rev. Lett. 78, 2275 (1997)], although Stratonovich [R. L. Stratonovich, Izv. Vyssh. 
Uchebn. Zaved., Radiofiz. 8, 116 (1965); Probl. Inf. Transm. 2, 35 (1966)] considered this quantity 
already in the mid-1960s. 
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under an arbitrary rotation in this subspace. (iii) As a result the reduced density matrix 
pi is also invariant under rotation in the subspace of subsystem 1 with S\ = N\/2, and is 
proportional to the identity matrix in this subspace. Thus 



A' l 



J- 1 



Ml 



2 



and Ei = log(iVi + 1) (in agreement with Ref. [32]). Similarly E2 = log(iV2 + 1)- Thus 

E = [log(iVi + 1) + log(iV 2 + 1) - \og(N + l)]/2. (32) 

We find in both the equal partition (N± = N2 = N/2) and unequal partition (N\ <C ^2 = 
N — 7V~i) limits, the entropy grows logarithmically with subsystem size iVx, 

lim E^ilog^). (33) 

yVi^oo 2 

3. ENTANGLEMENT ENTROPY OF ONE DIMENSIONAL SPINLESS FREE 
FERMIONS 

3.1. Fermionic Wave Function Approach 

3.1.1. Reduced Density Matrix of Spinless Free Fermion on a ID Lattice 

Now let us consider the application of this general formalism to the case of free fermions 
on a one dimensional lattice of infinite length. The subsystem is taken to N consecutive 
lattice sites. The Hamiltonian is simply: 

H = -t ^2(c j n+l c n + cic n+1 ). (34) 

n 

The Hamiltonian is simply diagonalized by performing Fourier transformation, let 



dke- ikn c k (35) 

\= E e ikn c n . (36) 

(37) 



'2ix 

n=—oo 
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We get: 



H = -2t dk cos kc\c k . (38) 



The ground state we will be interested in is to fill the vacuum state up to the Fermi 
momentum: 



i*o> = n 410) = (n(iE e "^)) i°) (39) 

k=—kf \ k n / 

But the mode of interest in our problem is the real space state, and according to our 
general formalism, we shall try to find out our 'state annihilation operator'. For a finite 
lattice of length N, any state can be written in such a general form: 

iK})=n^) n! i o )= a U}i°)' (40) 

i 

where i = 1,2, N, = {0, 1}, {n{\ denoting a set of occupation numbers. 

However, since the ground state \^q) is n °t vacuum state for q, the 'state annihilation 
operators' can be spanned as a s = YliOi, with Oi G {cj,c|}. Noticing that Wick's theorem 
is valid in this system, one can verify our assertion at the end of last section, and represent 
p s as 

Ps = tr e (p ) = ((EL O i ) t ) J] Oi, with 6 % e {ci, c\, c\c u c4} (41) 

{Oi} i 

Suppose we can find a unitary transformation: d m = Y2n=i v mnC n , so that (didj) = 
0, (fyo\d\dj\fy ) = 5ij{^o\dldi\^o). We can write po in terms of {di,d\}. Then according to 
Wick's theorem, and the fact (cjCj) = (d n d m ) = 0, only terms of G {d n d\, d\d n } survive. 
Then it is not hard to see that p s can be represented as a product: 

L 

P* = II ((Voldidll^dJl + (^o\dld0o)d\d^ (42) 

8=1 

3.1.2. Closed Form for the Entanglement Entropy 
According to equation 42, we see 
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i 

Follow the work of Jin and Korepin[24], we introduce Majorana operators of qs' 



C%-\ = Ci + C\ 



C21 = -i(ci - c\) 



dke~ lk \c k + c]' 



dke 



-ikn-(~\ 



Ck) 



dii-i — di + d\ and d 2 i = —i(di - d\). 



They satisfy: 

q = c], {ci,c m } = 26i m , (ci) = 0. 
The two-point correlation functions now are: 



\C21-1C2m-1) 



(c2l-lC2m) 



So we have: 



1 

2^ 
1 

2^ 
1 

2^ 
—% 
2^ 



dhdk^+^i^+ijih.+ci)) 

dke -id-m)k = = {c2iC2m) 

dhdk 2 e- ik ^ m ((c kl + cl)(-i)(c k2 - 4 2 )> 

dke -iW-m) g{k) = -fac^^) 



1 if k > kp or k < —kp, 
— 1 if kp > k > —k F 



here B L = 



(ciCm) = Sl m + i(B L ) lm , 



1 
-1 



with Gl 



go g-i ••• gi-L^ 
91 go i 



\9l-i go J 
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where g\ is defined as: 



.1 if k > kp or k < —kp, 
9i = — I dke~ Uh g(k), g(k) = { (48) 

-1 if kp > k > —kp 



2iv 



As we found, once we diagonalize the correlation-function matrix, the entanglement en- 
tropy is then given by — Vi lni^, where Vi is the ith eigenvalue of the matrix. However, 
due to the convention we adopted for the Majorana operators, the transformation is not 
unitary, but with an extra factor of 2. Thus, in order to ensure the normalization condi- 
tion of the density matrix, we must divide this factor of 2. Then the closed form for the 
entanglement entropy is readily given as: 

L 

s{ Ps ) = Y,< 1 ^) ( 49 ) 

i=l 

with 

e(x, u) = — — ^ — ^ n (~~^ — ) — ~~~2 — ^ n ^~~2 — ^' ^ S are e ^S envames °f Gi,. (50) 
However, to obtain all eigenvalues of v m directly from the matrix Gl is nontrivial task. 
Let us introduce 

D L {\) = det(Gx(A)) = det(AJ L - G L ) (51) 

II is a identity matrix of dimension L. It is known that Gp is a Toeplitz matrix (see [8]), 
i.e. its matrix elements depend solely on the difference between the two indices. Obviously 
we also have 

L 

D L (X)= H(\-v m ). (52) 

m=l 

From the Cauchy residue theorem and the analytical property of e(x, v), S(p s ) can be 
rewritten as 

S(p a ) = lim lim — / e(l + e)d\nD L (\). (53) 
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R = 8/2 




X = x+'\y 




FIG. 2: The contour c(e, <5). Bold lines (— oo,— 1 — e) and (1 + e, oo) are the cuts of integrand 
e(l + e, A). Zeros of Dl(\) [51] are located on bold line (-1,1) and this line becomes the cut of 
dlog Dl(X) for L — > oo. The arrow is the direction of the route of the integral we take and R is 
the radius of the circle. 

Here the contour c(e, 5) in Fig.l encircles all zeros of D L (X), but the function e(l + e, A) 
is analytic within the contour. The Toeplitz matrix Gx(A) is generated by the function g(k) 
defined by 



(54) 



A — 1 —hp < k < k F 

A + 1 k F < k < {2ix - k F ) 

Once the determinant of the Toeplitz matrix Gl is obtained analytically, one will be able 
to get a closed analytic result for S(p s ). 



3.1.3. The Toeplitz Matrix and the Fisher- Hartwig Conjecture 



Toeplitz matrix XL[0] is said to be generated by function (f)(9) if 



TM = (0. 



1-3 J i 



i,j = 1, ■■ ■ ,L - 1 



(55) 



where 



2tt 



2tt 



e~ iW d9 



(56) 
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is the l-th Fourier coefficient of generating function (f)(9). The determinant of Tl[0] is denoted 
by D L . 

Fisher-Hartwig Conjecture: Suppose the generating function of Toeplitz matrix (f)(6) is 
singular in the following form 

R 

(f)(9) = m\{tp T ,e r (e)ua r> eM (57) 

r=l 

where 

tpr,eM = exp[-^ r (7r -9 + 9 r )}, 9 r <9<2n + 9 r (58) 
u arj e r (9) = (2 - 2cos(# - 9 r )y\ Ua r > ~ (59) 

and if): T — > C is a smooth non-vanishing function with zero winding number. Then as 
n — > 00, the determinant of Tl[0] 



Dl = (m) L (jl^^J £ W> fo}> {A}, {9i}\, L - 00. 



(60) 



Here J-'l'ip] = exp f o 27r ln^(0)d0). Further assuming that there exists Weiner-Hopf fac- 
torization 

i)(9) = Tty] i, + (exp(ifl)) i). (exp(-ifl)) , (61) 
then constant S[if>, {«;}, {Pi}, {9i}} in Eq. 60 can be written as 

R 

S[if),{ai},{Pi},{9i}} = £[4>]l[G(l + a l + p i )G(l + ai -p i )/G(l + 2a l ) 



i=l 

X 

i=l 



J] (exp(i^)) J [*!>+ (exp(-i^)) J 
n l-exp(i(^-^)j , (62) 



\<i^j<R 



G is the Barnes G-function, £[ip] = exp(J^^l 1 kskS-k), and is the fc-th Fourier coefficient 
of In if) (9). The Barnes G-function is defined as 



00 

,2; 



G(l + z) = ( 27 r) 2 /V (2+1)2/2 -^ 22/2 Y[{(1 + z/n) n e- z+z2/ ^}, (63) 

71=1 

where 7^ is Euler constant and its numerical value is 0.5772156649 • • • . This conjecture has 
not been proven for general case. However, there are various special cases for which the 
conjecture was proven. 
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For our case, the generating function g{9) has two jumps at 9 = ±kp and it has the 
following canonical factorization 

m=mtfhw,k F mfhw,-kAo) ( 64 ) 

with 

lW) = (A + 1) (j^yj , 0(A) = -fr(A) = &(A) = — ln^. (65) 

The function t was defined in Eq. 58. We fix the branch of the logarithm in the following 

way 

-7T < arg (^jy) < n - ( 66 ) 

Here we verify the the above factorization expression explicitly. We have two jumps at 
9 = ±k F . According to Eq.33, we have two t function: 



taukw = expf— i/3i(7r - 9 + kp)], kp < 9 < 2n + k F 

(67) 

tf3 2 -k F = exp[— 2/3 2 (7r — 9 — k F )], —k F < 9 < 2n — hp 



Then g{9) is given as: 



9(0) — (A + 1) ( T - — T ) expH(l) In(^)(7r - 6> - Ml 







~k F / 


(r 






-1 




+ 1 
- 1 






-k F / 


(r 







-9 + k F )] (68) 
n — 9 — kp A + l., r 7r — 9 + kp A + l., 

< A+1 > li^rj exp| — — 1ii( a3t)] ^-^r~ ,n{ —i )] 

When 6> G [—k F , k F ], 

A + 1 n-e-k F 

. However, for t$]. F , 9 is not defined in this region, but considering periodicity, 9 G [2n — 
kp, 2n + kp], and to move it into the same region as of tp 2 _ kp , we must write tp u k F as 

, A + 1 > K-(9+2ir)+k F A + 1 kp-B-Tc 



t/h,k F 



2tt 



A-r v A-r 

Then we have 
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g(0) = (A + ^(^Ll)" 1 = A - 1, when 9 e [-k F , k F ] (69) 
A — 1 

Similar argument also applies to 9 G [kF, 2vr — hp}- And we have verified that the factor- 
ization works. 

For A i [-1,1], we know that |s»(/?i(A))| < \ and |9£(/? 2 (A))| < | and Fisher-Hartwig 
conjecture was PROVEN by E.L. Basor for this case [4]. Therefore, we will call it the 
theorem instead of conjecture for our application. Hence following the theorem in Eq. 60, 
the determinant -Dl(A) of AIl — Gl can be asymptotically represented as 

-/3 2 (A) 



D L (X) = (2-2cos(2A; F )) {c(l + /3(A)) G(l - /3(A)) } 

(A + l)((A+l)/(A-l))^ F/ i L L- 2 ^). (70) 



Here L is the length of sub-system A and G is the Barnes G-function and 

G{1 + /3(A))G(1 - /3(A)) = e~( 1+ ^ 2 W f[ { (l - Vw/» a j . (71) 



Therefore, 



lnDi(A) = L(ln(A + 1) - — ln(A + 1) + — ln(A - 1)) - 2/3 2 (A) In I 

7T 7T 

+ 2(ln(G(l + (3)) + ln(G(l - /3))) - /3 2 (A) ln(2 - 2 cos 2A; F ) 



^ L( ^ + ^_^ w(A)inL 
- 2 <ra + ^™'"< 2 -« 



(72) 



(73) 



Here our only concern is the first term and the second term which diverge linearly or 
logarithmically as the size the subsystem. 



3.1.4- Asymptotic Behavior of the Entanglement Entropy 

Now let us proceed to calculate the leading order of the entanglement entropy according 
our results. 
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First, we consider the term grows as subsystem size L. It is not difficult to see that the 
contribution from this term is actually zero, since the only residues arising from poles at 
A = ±1 are just zero. 



S lM = lim Urn — La P +C + W +t + A ) - 1 + £ - A H 1 + C ~ X )) 
1KHsJ e ^ + 5^0+ 2?rz J c V 2 2 2 2 J 

1 -fc F /7T _ fc F /7T (74) 

1 1 + A 1-X 1 
= 

Now let us turn to the second leading term of order In L 



2 r e(lHAMA) 
52 = ^/ dA (l + A)(l-A) lnL 

2 /* dA / 1 + 6 + A, ,1 + e + A, 1 + e-A, ,1 + e 



r ^ ( — ln( - ) - — ln(— ) ) (75) 



1 , + InL 
x In 



2m \X-lJ (1 + A)(1-A) 
This contour integral can be calculated as follows. First, let us look at Fig.l, noting that 



I dA(...)= ( [+ f+ f+ MdA(...), 

Jc(e,8) \J AF JFED J DC JCBA/ 

according to contour integral theory, this contour yields the same result as 



(76) 



dA(...)= (/ + / + / + / ) + (/ + / ) ) dA(...). (77) 

c(e,5) V J AF JFD J DC JCA J FEDF JACBA / 

The contour j FED ^ and j ACBA ' are merely closed circles around the points ±1. As we 
take the limit e — > + , they shrink and contain only this two points. Using Cauchy's residue 
theorem, it is easy to see that the contribution from these two contours is zero. 

So the contour integral is simplified to 

* = 5? (i. 1+ ,„ t - 7_ 1+ „-J dA (l + A)(l-A) < 78) 
The rest part of the integral function is analytic within the contour, however, /9(A) could 
have jumps in its angular part since we fix the branch by requiring —n < arg (in jzj)- 
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8(x + zO ± ) = In 

= 2^( v ln T^ +ln rrri^J 

--L^ilf+lnHl + tO^))) 

\ 1 — X / 

= 2S ('"S + ln(e ' ( ' +0±)) ) ' 
Taking into account the branch-cut condition, one immediately see that 

l ne *(T+o + ) = _ z(7r _ +) 

l ne '(T+o-) = 1(^-0") 

Therefore, our contour integral can now be written as 

2 /•' e(l,A) , 1 

S2 = ^1,^(17100^) 5 

Thus we have obtained the leading order behavior of the entanglement entropy of a 
segment of length L embedded in an infinite spin-less free fermion lattice. Here we did not 
calculate the sub-leading terms, however, when hp, the Fermi momentum or filling factor, 
becomes extremely small or very close to 1, they will become more and more important 
and eventually kill the entanglement's logarithm dependence on L. The criterion is given by 
C = 2L sin kp > I. 

3.2. Bosonic Approach 

3.2.1. Brief Introduction to Bosonization of ID Spinless Fermion 

In this section, we shall briefly introduce the bosonization of fermionic systems in one 
dimension. We will generally follow [38] which follows Haldane's constructive approach. 

For simplicity we will only consider the bosonization of a theory involving only one species 
of fermions. And bosonization of a theory is possible whenever the following prerequisites 
are met: 
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(79) 



(80) 



1. The theory can be formulated in terms of a set of fermion creation and annihilation 
operators with canonical anti-commutation relations: {c^cj,,} = 5kk', k G [—00,00]; 

2. The label k above is a discrete, unbounded momentum index of the form k = -^-{n^ — 

with n k G Z and G [0,2). Here n k are integers, L is the length of the system 
size, and 5b is a parameter that will determine the boundary condition for the fermion 
fields defined below. 

Obviously the prerequisites here are not directly satisfied by the lattice free fermion model 
we considered in previous sections. However, they can be satisfied by doing the following 
procedures: i)performing a particle-hole transformation c\ — > for k < —hp, ii) then 
shifting the Fermi point to kp = 0, iii) letting the lattice spacing a — ► 0, iv) extending 
definition of k to (—00, 00). 
The fermion fields: 

= (f) 1/2 E ( 82 ) 
fc=— 00 

with inverse c fc = (2vrL)" 1/2 / dxe^^x). (83) 

J-L/2 

And given a set of discrete fc's of the form above, the fields ip(x) obey the following 
periodicity condition: 

tp(x + L/2) = e inSb tp(x - L/2). (84) 
Using the following identity [18] 

J2e iny = 27r^%-27m) (85) 
we can immediately get the anti-commutation relations: 



{^(x), V>V)} = Y E e- i(x - x ' )(n - <5b/2)2 " /L = 2tt ^ <5(x - x' - nL)e innSb ; (86) 
{V(aO,^(aO} = °- (87) 
Vacuum State |0) : 
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Let |0) be the state defined by the properties 



c fc |0) = for k > 0, (88) 
4|0) = for k < 0. (89) 

We shall call |0) the vacuum state and use it as reference state relative to which the oc- 
cupations of all other states in Fock space are specified. With this definition we can define 
the operation of fermion-normal-ordering, to be denoted by * *, with respect to this vacuum 
state: to fermion-normal-order a function of c and c^'s, all c k with k > and all c k with 
k < are to be moved to the right of all other operators(i.e. all all c\ with k > and all c k 
with k < 0), so that: 

IABC ...: = ABC ■ • • -o (0\ABC . . . |0) for A, B, C, ■ ■ ■ e {c k] c\} (90) 

N -particle ground state \N) Q : Let N be the operator that counts the number of electrons 
relative to |0) : 

oo oo 

N= * c l c **= E [ck-o(0|4c fc |0) ] (91) 

k=—oo k=— oo 

The set of all states with the same iV-eigenvalues N will be called the iV-particle Hilbert 
space Hn- It contains infinite number of states, corresponding to different particle- hole 
excitations. Let us denote all of them by \N). For a given N, there is a state which contains 
no particle-hole excitations. We will denote it as \N) Q . To prevent possible ambiguities in 
its phase, we define it by the following order: 



Bosonic operators b q and b q : 



4c[ JV _ 1) ...c}|0) foriV>0, 

|0) foriV = (92) 

^c N+1 c {N+2 ) • • • c |0) for N <0. 



oo . oo 

I s. — > + . —I 



v ^ k=— oo v k=— oo 

with q = ^n q > where n q G Z + is a positive integer. Thus the bosonic creation and 
annihilation operators are defined for q > only. 
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And it is not hard to prove the following bosonic commutation relations: 



K M = K b U = °> Wf> V] = [N 9 , bl = 0], for all q, q'; (94) 

OO j 

[ 6 9> 6 J'] = — ( C l+ g - 9 ' C fe -4+^+9') = V- ( 95 ) 



n 9 

fc = — OO ^ 



Making a connection with Eq.(93), we could see that in each iV-particle Hilbert space 
Hn, \N) functions as vacuum state for bosonic operators defined above: 

6 g |iV) = 0, for all k. (96) 

This is because \N) is the iV-particle ground state and does not contain any particle-hole 
excitations(bosonic excitations) . 

With a proper construction of a set of bosonic operators in k space, the construction of 
boson fields is just straightforward: 

- —e- igx b q e- aq/ \ tp\x) = ~J2 —e iqx b\e~ aq/2 . (97) 



<p{x) 

Here a > is an infinitesimal regularization parameter which is used to regularize q — > 
oo divergent momentum sums that arise in certain non-normal-ordered expressions and 
commutators. The following commutation relations can be verified for the boson fields 
defined above: 



[<p{x),<p(x')] = i^(x),^(x')]=0, (98) 

[ip{x), ip ] (x')} = — e - q[i{x - x ' )+a] = - ln[l - e -<£(*-*'-<«0] (99) 
q >o Uq 

Then consider their Hermitian combination: 

(j){x) = <f(x) + <f t (x). (100) 
Check the canonical commutation relation: 



{x),d x «f)(x')] = i^- ( -2- — ]— h 



Xj \ i-£{x-x'-ia) g— i*j-{x— x'+ia) ^ 



(101) 
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which is exactly the canonical commutation relation for boson fields when L — ► oo. 
Bosonization and entanglement entropy 

Let us look at the corresponding mapping in real space since the partition of the system 
is usually carried out in real space. Therefore, to legitimate our use of bosonization to study 
entanglement entropy, we have to establish the relation between the fermion fields and the 
boson fields in real space and show that the mapping, at least approximately, preserve the 
partition of the system. 

First, let us look at the fields: 



^ = E ^ e ^ e ~ a<z/2 ^= E c U c * 

q>0 V 1 V 1 k=-oo 

-? ^ i r L/2 r L/2 

T^ n q , 27tL J -lit. J-L/2 



dx x dx 2 ^(x 1 )ip(x 2 ) e l{ - aq/2 - qx+qXl - kxi+kX2) (102) 

q >o 27rLnq k 

dx 1 j^^ix^ixJe^^-^V 

n q >0 q 

dx 1 y^(x 1 )xlj(x 1 ) ln(l - e i2 -r {xi ~ x - a/2) ). 

Lj 

It seems that the mapping of fields does not fulfill our requirement. However, if we look at 
the fermion density operator which is the object one will directly work with when calculating 
entanglement entropy, we shall see that our requirement is indeed satisfied. 



2-7T 

Pf (x) = :^{x)^{x): = T Y1 e ~ iqx E * c U c ** 

q x 

J2i^n q (e- iq % - e iq %) + — "£>W* ( 103 ) 



L ^ v qy q L 

q>0 k 

= d x <p{x) + ^N, for a -> 0. 

Lj 

For the whole system, the particle number is conserved, so N is just a number. Thus we 
have justified our utilization of bosonization to study the entanglement entropy in many- 
fermion systems. 

However, we still need to find out corresponding bosonic states of the fermionic states we 
are interested in. At present we are only interested in fermion ground state at zero temper- 

25 



ature, i.e. the Fermi sea. It is the |0) we defined in Eq.(88) and (89) which corresponds to 
the vacuum state \N = 0) of the boson modes. 

3.2.2. Entanglement Entropy of free Bosons 

In this part, we shall first introduce available analytic result in lattice model. However, 
though this approach can give a nice analytic expression for entanglement entropy, we can not 
get desired result explicitly due to technique difficulty. Then we will give a short introduction 
of the field theory approach, following Calabrese and Cardy[9]. 

Lattice Model 

Let us consider a system of coupled harmonic oscillators in which the Hamiltonian can 
generally be written in matrix form as: 

H = y][TjjPiPj + VijXiXj], (104) 

where the operators Xj's and pj's obey the canonical commutation relation: = iSij. 

Consider its correlation functions of positions and momenta 

-^■n,m y^n-^m/ i Pn,m (PnPm) ■ (^0^) 

The ground state is a Gaussian state, and the multi-point correlation functions observe 
Wick's theorem: 

) = (x n x m )(x k xi) + (x n x k )(x m xi) + (x n xi)(x m x k ). (106) 

This also holds inside the subsystem when we do the truncation. According to Wick's 
theorem, this indicates the density matrix of the subsystem, i.e. the reduced density matrix 
is an exponential of momenta and spatial coordinates. Therefore, in principle we should be 
able to write the reduced density matrix as: 

p s = K,< "' = e -^MjPM+ v i^\ (107) 

here JC is a normalization factor. It is in general not easy to obtain an explicit analytic 
expression for X" and V due to two reasons: first, even in the simplest case of nearest 
neighbor coupling, terms that do not conserve particle numbers would arise; second, the 
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transformation must be symplectic 3 . However, for such a Hamiltonian, it is always possible 
to find a symplectic transformation Cj = -75(2^ +ipj), «fc = ^2j(Ak(j) c j + ^lU) c ] "which can 
symplecticly diagonalize the Hamiltonian to the following form: 

H' = J2 £ ( k ) a W- ( 108 ) 

k 

The eigenvalues e(fc)'s follow from the eigenvalues v\ of the X'P' A matrix via 

cth(e(fc)/2) = v k /2. (109) 
And the entanglement entropy is given by 



S = J2({"k + 1/2) ln(z/ fc + 1/2) - K - 1/2) ln(i/ fc - 1/2)). (110) 

k 

However, to obtain the full spectrum of the matrix X'P' and calculate the entanglement 
entropy is highly non-trivial. No pure analytic approach has been developed, but with the 
help of numerical methods, one still could reproduce desired scaling law of the entanglement 
entropy. In the strong coupling limit which corresponds to a massless free field, the scaling 
behavior has been shown indeed to be | log L. 
Field Theory Approach[9] 

Consider a quantum field theory in one dimension space and one time dimension, de- 
scribed by the following action 

S = J \ (-(tt) 2 + mV) *r = j\ {{d^f + mV) d\ (111) 

where ir(x) = id^ip is the canonical momentum, and satisfies the canonical commutation 
relation [ip(x),7r(x')] = i5(x — x'). The density matrix p in a thermal state at inverse 
temperature (3 is 

P{W)"}\{^)'}) = Z(f3r\{<f>(xr}\e- H \{<f>(x'y}), (112) 

3 "simplectic" here means preserving the commutator [xi,pj] = <5y, i.e. a transformation that preserves the 
" symplectic" matrix ^ 1 J ' 

4 Here ' indicates that these are obtained by truncating the original matrices X and P. 
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where Z(/3) = tre~^ H is the partition function, and {(f>(x)} are the corresponding eigenstates 
of ip(x): ip(x)\{(p(x')}) = 4>(x')\{(j)(x')}) . This can be expressed as a (Euclidean) path 
integral: 

p = Z- 1 J [dcj>{x, r)] J] 5{<j>(x, 0) - cj>{x')') H 8{cj>{x, (3) - cj>{x")")e- s \ (113) 

X X 

where Se — f$ LecLt, with Le the euclidean Lagrangian. The normalization factor Z, i.e. 
the partition function is found by setting {tp(x)"} = {ip(x)'} and integrating over these 
variables. This has the effect of sewing together the edges along r = and t = j3 to form a 
cylinder of circumference (3. 

The reduced density matrix of an interval A = (x{,Xf) can be obtained by sewing together 
only those points which are not in the interval A. This has the effect of leaving an open 
cut along the line r = 0. For the calculation of entanglement entropy, we need to perform 
a replica trick here. Instead of calculating trp^logp^, we compute p n A first, for any positive 
integer n. To do this, we make n copies of above set-up labeled by an integer k with 1 < k < 
n, and sew them together cyclically along the open cut so that ip(x)' k = ^(x)' fc / +1 (and ip{x)' n = 
i/j(x)i) for all x € A. Let us denote the path integral on this n-sheeted structure by Z n (A). 
Then 

^ = (H4) 



and 



<? r d Mil 

SU = -hm- — — . (115) 

n-*i dn Z n 



For a free theory on such n-sheeted geometry, it is easier to use the identity 5 : 

^log^ = -^ J G n (r,r)d\ (116) 

where G n (r, r') is the Green's function in the n-sheeted geometry. To obtain , we need 
G n — nG\. The Green's function can be obtained by solving the Helmholtz equation 



V r 2 + m 2 )G n (r,r') = 5 2 (r-r') (117) 



This only holds for non-interacting theories. 
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with n-sheeted geometry. In polar coordinates(2D) this simply means extend the domain of 
9 to [0, 27m). However, an explicit formula is only available for infinite volume, which means 
we have to extend the subsystem to a semi-infinite one. 

After solving for G n (r, r'), let r' = r and perform the integral, we have 

^log(§) = -\J (G.(r,r) - n Gl (r,r)) = ^(n - I). (118) 
This will lead us to the final entanglement entropy 

S = -trplogp = - Atrp»| n=1 = -A( m 2 a 2)^(n-l/n)| n=i = _J_ logm 2^ (ng) 

The power of a is inserted to make the result dimensionless, following Cardy's convention. 
In the massless limit where m — > which we are more interested in, physically it is natural 
to replace m with the inverse subsystem size L _1 since the coherent length diverges and 
subsystem size is the only relevant characteristic size in this problem. Also be careful that 
we actually dealt with a semi-infinite subsystem which has only one boundary of partition in 
above approach 6 . When we consider a more common subsystem which has two boundaries 
of partition, we want to double this entanglement entropy as it is a boundary effect. Thus 
we recover the result S — 3 log L in agreement with other approaches. 



4. NEXT STEP: GENERALIZATION TO HIGHER DIMENSIONS 

4.1. Violation of Area Law in Fermionic Systems in Higher Dimensions: Known 
Results 

In Wolf's work [40], he considers a general number preserving quadratic Hamiltonian 

H = Y J T a ,^ a c (3 ,T = T\ (120) 

a,f3 

describing Fermions on a d-dimensional cubic lattice, so that each component of the vector 
a, (3 corresponds to one spatial dimension. Peschel et al.[ll, 31] obtain a general result on 
the reduced density matrix of such a Hamiltonian 

p s = JCe~ H ' = Kj >-;•."',';<.. (121) 

6 This is verifiable in several cases. 
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where H' = log(l — C)/C, C = (c\cj)ij. Here all the indices could be vectors when we 
consider dimension d > 2. Their argument is quite general and applies to not only arbitrary 
dimensions but finite temperature as well. And from this p s one can easily derive the 
following expression of entanglement entropy: 

L d 

S = 52h(\j) t (122) 

j'=i 

,.. 1 + X , 1 + 2 1— X, 1 — X 

h(x) = ^ X °Z— 2 2~ l0g ^~' (123) 

which are general and also hold for systems of higher dimensions. Ay's are eigenvalues of 
the correlation matrix C. However, a direct computation of S via the diagonalization of C 
is highly non-trivial even in the simplest case as we have seen in previous sections. Wolf's 
argument is based on the upper bound and lower bound behaviors of the entropy function 
h(X). 

His conclusion is that 

c-L d - 1 log L < S < c+L A - x (log Lf, (124) 

with constants c± depending only on the Fermi sea. This result requires that the Fermi 
surface must be regular enough, i.e. not fractal nor Cantor-like. 

Gioev and Klich provide a more specific expression for the scaling behavior of fermion 
entanglement entropy by making a connection with Widom conjecture 

S ~ L (o u-i^ Tt I I \n x ■ n p \dS x dS p , . (125) 

[Z7T) LZ J an Jg r 

This has been verified numerically [2, 28]. 



4.2. Outlook: Bosonization in Higher Dimensions and Entanglement Entropy 

High Dimension Bosonization We have shown basics of one-dimensional bosonization in 
previous sections. Bosonization in arbitrary dimensions was first formulated by Haldane[20] 
(for a review see [22]). The basic idea of bosonization in d > 1 dimensions is to divide the 
Fermi surface into small segments S with height A in the radial direction and area A d ~ x 
along the Fermi surface. These two scale must satisfy the following condition: 

A; F >A>A. (126) 
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Then we focus on the low energy physics relative to Fermi energy. This is done by integrat- 
ing out the high momentum (energy) degrees of freedom to get the effective Hamiltonian. 
Working with the effective Hamiltonian within each segment S, we will find ourselves in 
a situation similar to that near the Fermi points in one-dimensional systems. For small 
momentum transfer q <C A, it is again possible to pair up the particle- hole excitations in a 
bosonic way as in one-dimension case within small correction as long as the prerequisites 
are fulfilled. 

Possibility of Entanglement Entropy via Bosonization 

The possibility of application of bosonization in higher dimensions arises from several 
aspects: 

1. According to results in one-dimensional system (given by Korepin et al.) where 
correction with respect to the Fermi energy Ef is included, we see that the singular behavior, 
i.e. the logL scaling, only emerges when Ef becomes big enough or the Fermi sea becomes 
deep enough. This is one of the most fundamental requirements for the 1-d bosonization to 
work properly. In one dimension bosonization becomes exact when we have a infinitely deep 
Fermi sea; 

2. The agreement on the scaling behavior of entanglement entropy in free fermion systems 
and free boson systems; 

3. There has been successful application of bosonization in the calculation of entan- 
glement entropy of two dimensional free Dirac fields[10]. Even though in that case the 
Fermi surface is absent, it is still nonetheless a strong indication that bosonization in higher 
dimensions could work in the presence of Fermi surface. 

The advantage of bosonization approach is that it could take into account the interactions 
in arbitrary dimensions. 
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